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Qh Abstract 

A class of explicitly integrable models of 1+1 dimensional dilaton gravity coupled to 



(N 



scalar fields is described in some detail. The equations of motion of these models reduce 
to systems of the Liouville equations endowed with energy and momentum constraints. 
The general solution of the equations and constraints in terms of chiral moduli fields 
is explicitly constructed and some extensions of the basic integrable model are briefly 
' discussed. These models may be related to high dimensional supergravity theories but 

here they are mostly considered independently of such interpretations. A brief review of 
' other integrable models of two - dimensional dilaton gravity is also given. 

o . 

1 Introduction 

It is well known that 1+1 dimensional dilaton gravity coupled to scalar matter fields is a 
O 1 reliable model for some aspects of high dimensional black holes, cosmological models and 

branes. The connection between high and low dimensions has been demonstrated in different 
contexts of gravity and string theory and in some cases allowed to find general solution or some 
special classes of solutions in high dimensional theories 1 . For example, spherically symmetric 
gravity coupled to Abelian gauge fields and massless scalar matter fields exactly reduces to a 
1+1 dimensional dilaton gravity and can be explicitly solved if the scalar fields are constants 
independent of coordinates. Such solutions may describe some interesting physical objects - 
spherical static black holes, simplest cosmologies, etc. However, when the scalar matter fields, 
which presumably play a significant cosmological role, are not constant, few exact analytical 
solutions of high dimensional theories are known. Correspondingly, the two - dimensional 
models of dilaton gravity nontrivially coupled to scalar matter are usually not integrable. 

To obtain integrable models of this sort one has to make some serious approximations, in 
other words, to deform the original two - dimensional model obtained by a direct dimensional 
reductions of realistic higher dimensional theories 2 . Nevertheless, these models may qualita- 

* Alexandre . FilippovQjinr . ru 
For a more detailed discussion of this connection and references see e.g. 0, |5] and next section 
2 Note that some important four - dimensional space - times having symmetries defined by two commuting 
Killing vectors may also be described by two - dimensional dilaton gravity. For example, the cylindrical 
gravitational waves are described by the 1+1 dimensional dilaton gravity coupled to one scalar field. The 
stationary axially symmetric pure gravity may be described by a 0+2 dimensional dilaton gravity coupled to 
one scalar field (this may be related to the previous cylindrical case by the analytic continuation of one space 
variable to imaginary values). Similar but more general dilaton gravity models were also obtained in string 
theory. Some of them may be solved by using modern mathematical methods developed in the soliton theory. 
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tively describe some special physically interesting solutions of higher dimensional gravity or 
supergravity theories related to the low energy limit of superstring theories. 

When preparing and writing this report I often recalled very stimulating talks with Misha 
Saveliev about the Liouville and other integrable theories. I cannot express my gratitude to 
him in this world but I wish to mention that he strongly influenced my interest and my 
approach to matters discussed below. 

2 Some integrable models of 1+1 dimensional dilaton gravity 
coupled to scalar matter 

The effective Lag rangian of the X — I — X dimensional dilaton gravity coupled to scalar fields ip n 
that may be obtained by dimensional reductions of a higher dimensional spherically symmetric 
(super)gravity can usually be (locally) transformed to the following form 3 



£ (2) 



ipR(g) + V(<p, il>) + %n{v, VO 9 ij diip n d^ n 



(1) 



Here gij(x°,x 1 ) is a generic 1+1 metric with signature (-1,1), g = det|gjj| and R is the Ricci 
curvature of the two - dimensional space - time, 

ds 2 = g lJ dx i dx 3 , = 0,1). (2) 

The effective potentials V and Z n depend on the dilaton tp(x , x 1 ) and on N — 2 scalar fields 
t/j n (x°, x 1 ) (note that Z n < 0). They may depend on other parameters characterizing the 
parent higher dimensional theory (e.g. on charges introduced by solving the equations for the 
Abelian fields). Here we mainly consider the simplest kinetic terms, with Z n (ip,ip) = Z n ((p), 
or even with constant Z n that are independent of the fields. We also used in JI) a Weyl 
transformation to exclude the gradient term for the dilaton. 

To simplify derivations we will use the equations of motion in the light - cone metric, 
ds 2 = —4:f(u,v)dudv. By first varying this Lagrangian in generic coordinates and then 
going to the light - cone ones we obtain the equations of motion 

d u d v <p + fV(<p,iP) = 0, (3) 

fdi{di<p/f) = £ Znidi^n? , (i = U,v) . (4) 

d v (Z n d u i/j n ) + d u (Z n d v if) n ) + fV$ n (tp, ip) = Y z m^ n d u ip m d v ip m , (5) 

771 

d u d v In |/| + fV<p{tp, tp) = Y Z n,ip duipn d v ip n , (6) 

where V v = d v V, V^ n = d^ n V, Z n>ip = d v Z n , and Z m ^ n = d^ n Z m . These equations are not 
independent. Actually, © follows from (jSJ) — ©. Alternatively, if ©, Q, ® are satisfied, 
one of the equations (JSJ) is also satisfied. 

Consider the solution with constant scalar field ifj = ■i/'o (the 'scalar vacuum' solution). 
This solution exists if V^(ip,tpo) = 0, see eq.©. The constraints @ can now be solved 
because their right-hand sides are identically zero. It is a simple exercise to prove that there 

3 For a detailed motivation and specific examples see 0, where references to other related papers may be 
found. Due to space limitations, only absolutely necessary references are given here. 
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exist chiral fields a(u) and b(v) such that ip(u, v) = <p(r) and f(u, v) = <p>'{r) a'(u) b'(v), where 
r = a(u) + b(v) (the primes denote derivatives with respect to the corresponding argument). 
Using this result it is easy to prove that © has the integral ip' + N(<p) = M, where N(tp) 
is defined by the equation N'(ip) = V((p,ipo) and M is the constant (integral) of motion. 
The horizon, defined as a zero of the metric h(r) = M — N(ip), exists because the equation 
M = N(<p) has at least one solution in some interval of values of M. These solutions are 
actually one - dimensional ('automatically' dimensionally reduced) and can be interpreted as 
black holes (Schwarzschild, Reissner Nordstr0m, etc.) or as cosmological models. Actually, 
the scalar vacuum theory is equivalent to the pure dilaton gravity 4 , which is a topological 
theory, and this explains this drastic simplification of the equations of motion. These facts 
are known for long time and were derived by many authors using different approaches (for 
references see e.g the recent review j2]). 

The properties of the general dilaton gravity theories (^Q) are much more complex. They 
do not reduce to one - dimensional models and are in general not integrable. If Z n ^ the 
theory may be integrable only with very special potentials V(y>,ip) and Z n (ip,ifi). Roughly 
speaking, the two - dimensional models may be integrable if either the potentials Z n can 
be transformed to constants or the potential V is zero. In the first case, V should have 
very special form that is described below. In the second case, for the theory to be explicitly 
analytically integrable 5 the potentials Z n must be very special functions. 

In search for integrable 1+1 dimensional models it may be useful to also consider 'mixed' 
cases. Roughly speaking, this means the following. Suppose we have three different species 
of scalar matter fields - t/j n , Xki an d cr s . The characteristic properties of these fields are: the 
potential V depends on <p and ip n ; Z n and Z^ are constants, say, Z n = Zj, = — 1; Z s = Z s {tp) 
do not depend on a and ip. Thus the Z-part of the Lagrangian (^Q) is 




. n k s 



Looking at the equations @ - © we see that Xk are free massless fields (i.e. d u d v \k = 0), 
the equations ®, © for ip, and ip n are independent of the other fields, and, finally, the 
equations for a s are linear. The right - hand sides of the constraints @, 

J2Zn(d i ^ n ) 2 + J2 z k(diXk) 2 + J2 Z ^)( d i C7 s) 2 , (i = u,v), (8) 
n k s 

depend on all three species of the matter fields, so that the fields are not decoupled (this 
is natural because some linear combinations of the constraints give the total energy and 
momentum of the system that are constrained to zero) . The most general system of this sort 
may be explicitly integrable if its nonlinear sector is explicitly integrable and the constraints 
can be explicitly solved with arbitrary solutions for a and x fields. Unfortunately, in such 
a general setting the constraint cannot be solved. As is shown in the next section, when 
we have only fields /, ip and ■0 (by taking constant a fields), each constraint involves only 
functions of one variable (n or v, resp.). However, this condition cannot be satisfied for 

4 If V = V(ip) and Z n = the theory Q is called the pure dilaton gravity. It is integrable with arbitrary 
potential V{ip). 

5 We call the theory explicitly integrable if one can explicitly write its general solution analytically, in terms 
of a sufficient number of arbitrary functions, in our case, of chiral free fields and of their derivatives. For 
further use in physics problems we even need to have simple enough expressions for the general solution. 
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the term Z s {ip) (d{ a s ) 2 in the integrable models discussed here. For this reason, mainly the 
models with the matter fields ip and x wm be considered here. Nevertheless, it is useful 
to keep ex-fields in mind for two reasons. First, in some cases these terms may be taken 
into account as a perturbation. The second reason is the following. In further reductions 
to dimensions 0+1 or 1+0 the a fields may be exactly solved giving additional terms to the 
potential V. The resulting theory may be integrable or approximately integrable. We do not 
consider such possibilities here and concentrate on the theories with the fields ip and x- 

Let us first discuss the second 'pure' case (including x-fields is trivial and we omit it). So, 
suppose that V = 0, ip = x = 0, and Z s = Z(tp). Then, from the first equation we have (up 
to the coordinate transformation u i— ► a(u), b t— > b(v)): (p = u + v = r and the scalar fields 
satisfy the linear equation (t = u — v) 

(9 r 2 -^(M) + ^t^(M)=0. (9) 
Z(r) 

When Z(<p) = —ip this is the Euler - Darboux equation, for the general solution of which one 
can write a rather complex integral representation (see e.g. |3|). Moreover, the expression 
for the metric is very difficult to analyze. In general this does not allow us to find an explicit 
analytic expression for the metric. Thus, such theories may be called integrable but I would 
say that they are not explicitly integrable. However, if Z n (ip) are solutions of the one - 
dimensional Liouville equation, the potentials produced by such Z n (r), are 'non reflecting', 
the equations for a are explicitly integrable and are simple functions of the chiral free massless 
fields a(u) and b(v) (see [I]). The metric is then explicitly expressed through these chiral fields 
and thus the theory is explicitly (and elementary) integrable. These solutions may describe 
plane waves of scalar matter coupled to gravity or certain inhomogeneous cosmologies. 

The qualitative difference between the solutions of © for 'non reflecting' potentials and 
the 'realistic' ones (corresponding to Z = —tp) can be made more clear if we use the following 
beautiful, compact solutions of eq.@ with Z = — ip 2X+1 (A > — 4). The first solution was 
essentially found by S.D.Poisson, 

a(r,t) = / da q(t — r cos a) sin 2A a , (10) 
Jo 

the second one was derived later by E.W.Hobson 6 , 

/•CO 

a(r,t) = / da q(t ± r cosh a) sinh 2A a , (11) 
Jo 

where q(x) is an arbitrary 'generating function'. Now, the potential corresponding to Z = 
—ip 2X+1 becomes 'non reflecting' if A = I + i, with integer /. Indeed, then Q is equivalent to 

(d 2 r - d 2 )a(r, t) - l ^±^a{r, t) = , (12) 

where a = r l+1 a. The potential 1(1 + l)/r 2 is a special case of the 'non reflecting potentials' 
considered in 4 , with which the solution of this wave equation can be explicitly written as a 
linear combination of a free massless field x = a (t + r ) + b(t — r) and of its derivatives (with 
coefficients depending on r). One can easily prove that the solutions (|10j) and (jlljl reduce to 

6 Known from 19-th century, these solutions and their generalizations may be found in the reprint of a still 
useful old book 0. 
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such local fields when A = I + | with integer I. If I is non integer the solution is a continual 
superposition of waves with different 'velocities' sec a or cosh -1 a. In this sense, the solution 
is generally nonlocal and, in fact, not simpler than one used in It may be instructive to 
rewrite Ea, ()l(J|) . with half - integer A, in the form 

1 r t+r , 
a (r,t) = - dxq(x)\(x - (t - r))((t + r) - x)] . (13) 
r Jt-r 

If I = we find that 

a(r,t) = -[Q(t + r)-Q(t-r)} , 
r 

where Q'(x) = q(x). Similarly, one can find the local expressions of a in terms of massless 
free fields for any integer I. 

Originally, the equations (|l(Jj) and (|llj) were applied to description of simple cylindrical 
waves (see e.g. (Hj). More recently, it has been shown that they can be generalized to describe 
rather complex cylindrical gravitational waves (see e.g. |S| in which a generalization of these 
classical solution in the context of the Penrose twistor program is also discussed). In fact, 
the simplest cylindrical gravitational waves derived by A.Einstein and N.Rosen (see [7]) may 
be found by solving just the equation Q. To find the most general cylindrical gravitational 
waves one has to solve the following nonlinear equation 

dtirM^dtM) - <9 r (rM -1 d r M) = , (14) 

where M is a 2 x 2 symmetric matrix depending on the cylindrical metric. This equation 
may be solved by a generalization of the Poisson formula Q or by more general methods 
developed in the theory of solitons. 

In fact, the equation (|14j) describing cylindrical gravitational waves is a special case of 
the field equations of cr-models coupled to gravity. This equation and its generalizations 
may be attempted to solve by the inverse scattering methods (ISM) of the soliton theory. 
Applications of these methods started with papers [E] and [0] which consider integrability 
of the stationary axial gravity. Further results and references to numerous papers dealing 
with this subject may be found in JHJ and The methods for solving the integrable two 

- dimensional theories reducible to solving equations like (j!4j) are not elementary and, for 
this reason, are not included in our classification of the matter terms (|15|) that may lead to 
elementary solvable two - dimensional dilaton gravity. In our notation, the models that may 
be integrated by ISM have V = 0, Z n = = but, in general, more complex structure of 
the d-terms: 

^Z a p(ip;a)di(T a djap (15) 

In simple cases one may transform these terms to the diagonal form similar to (|15|) but 
with the potentials Z depending on u-fields. For example, the Lagrangian (2.8) of the two 

- dimensional dilaton gravity derived in ^5] for the stationary, axially symmetric Einstein 
gravity can easily be transformed to the form of our 'general' dilaton gravity Q but this 
does not help to solve the equations of motion by elementary methods. On the other hand, 
we will see that there exist reasonably wide class of physically interesting, integrable two - 
dimensional theories with simple (constant!) Z n but rather complex potentials V. 

The best studied examples of such integrable models (belonging to the first 'pure' class, 
with constant Z n ) are Callan - Giddings - Harvey - Strominger (CGHS) model (V = go ) 
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and the Jackiw-Teitelboim (JT) model (V = gup). In CGHS model R = and in JT model 
R = —g\. A strong generalization of these models, in which the two - dimensional curvature 
is not conformal constant (V = g+e 9 ^ + g-e~ 9Lp ) was proposed by the present author in 
Like CGHS and JT models, this model belongs to the 'mixed' case having free fields Xk- 
It was called the 2-Liouville model because the fields rp± = In / db <p satisfy two Liouville 
equations. The Liouville equations are, of course, explicitly solved in terms of two pairs of 
the chiral free massless fields a±(u), b±(v). By using the constraints (0} one of the free scalar 
fields x can be expressed in terms of a, b and of other x-fields. So one may be inclined to 
decide that the two - Liouville model is explicitly solved in ^2]. However, there remained 
one 'small' problem - the right - hand sides of the constraints are negative definite while the 
left - hand sides may change sign. In [13 I did not find an explicit analytic expression for the 
chiral fields a and b that fully satisfy the constraints. The explicit solution of the constraints 
(flf) was constructed only five years later when I have found a much more general iV-Liouville 
model ^1], pp. Here I will show how this class of models may be constructed and give its 
complete solution in terms of chiral moduli fields. 

3 N-Liouville 1+1 dimensional model and its solution 

Let us suppose that the theory is defined by the Lagrangian, 
£ (2) = 



tpR(g) + V(ip, ip) +g lJ Z n di^) n djtp n + J2 Zk diXk djXk 

V n k / 



(16) 



with the following potentials: 

N 

\f\V=J2^9ne qn , Z n = Z k = -l. (17) 

n=l 

Here / is the light - cone metric, ds 2 = —Af(u,v) dudv, and 

N N 

q n = F + a n lf + Y i'mamn = X! ^mO-ran , (18) 

771=3 771=1 

where ipx + ^2 = hi|/| = F (f = ee F , e = ±1), ib\ — tp2 = <f and thus a\ n = 1 + a n , 
a2n = 1 — On- By varying the Lagrangian (|16|) in N — 2 scalar fields, dilaton, and in g$j and 
then passing to the light - cone metric we find N equations of motion for N functions ip n , 

N 

tnd u d v 4>n = ^9m,e qm a mn ; ej = -1, e n = +1, if n > 2 , (19) 

777=1 

as well as two constraints, 

N 

d = fdi(d^/f) + Y(^n) 2 = m(i), i = (u,v), (20) 

71=3 

where 4Dj(i) = — J2i®i Xk) 2 is the contribution of the free fields % (see eq.©) that depends 
on one variable because Xk = «fc(^) + /3fc(f) (here and below k = 1, ...,K). 
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With arbitrary coefficients a mn these equations of motion are not integrable. However, 
as proposed in [TI], the equations ("i""))) are integrable and the constraints (f*"U)l can be solved 
if the iV-component vectors v n = (a mn ) are pseudo - orthogonal. Then the equations ffTTTfl 
reduce to N independent, explicitly integrable Liouville equations for q n , 

d u d v q n - g n e q - = 0, (21) 

where g n = e\ n g n , X n = J2 fmflLi an d £ = //|/| (note that the equations for q n depend on 
e n and a mn only implicitly, through the normalization factor X n ) . 

The expression for the original fields in terms of the Liouville fields q n may be found by 
using the orthogonality relations for a mn (to 7^ n) and the definition of X n = (for to = n) 
combined in the equation: 

N 

eiai m ai n = X n 5 mn = Tn^mn • (22) 

Z=l 

This equation may be written in the matrix form, if we define the matrices a = (a mn ), 
7 = (7m <W) an d e = (e m 5 mn ): 

rp -1 rp 

a ea = 7 , aja = e. (23) 
It is important to note that in addition to these relation the matrix a satisfies two conditions 

«ln = 1 + CL n , CL2n = 1 — CL n . 

Now, using ((2*2*]) or (|2*3*)) we can invert the definition (|18|) and get 

N N N 

VVn = e m ^2 a mnlnq n , F = -2 ^2 j n a n q n , f = -2 ^2 JnQn ■ (24) 

n=l ra=l n=l 

Writing the first equation for n = 1,2, expressing a\ n , am in terms of a n and using linear 
independence of the functions ip n we may find very useful identities ('sum rules') for 7 n ,a n 
and a mn for to > 3. We have 

N N N N 

Ipl = ~1pl + °") 2 ~ ^2 T^C 1 ~ a n) - X! ^ m a mn(l + On)7n j ( 2 ^) 

n=l n=l m>3 n=l 

AT N N 

1p2= Ipl X! 7*0- ~ a n) + V>2 ~ a ™) 2 + Y ^ m Y a ™n{t ~ a n )7n ■ (26) 

n=l n=l m>3 n=l 

Thus we immediately find that the following identities should be satisfied 

N N I N N N 

Y = °' Y ^ nCLn = ~ 2' ^ nCt " = °' a m«7n = a "m a ™7n = 0, TO > 3 . (27) 

n=l n=l n=l n=l n=l 

Also, using the orthogonality relations (|22[) it can be proven that one and only one norm 
7 n is negative. We thus choose 71 < while other j n are positive. In physically motivated 
models the parameters a mn , j n and g n may satisfy some further relations. For example, the 

7 Here we suppose that An ^ and g„ 7^ 0. Otherwise the solution of the constraints should be modified in 
a fairly obvious way. We also denote j n = A n ~ ■ 
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signs of 7 n and g n may be correlated so that g n /ln < 0. However, such relations do not 
follow from the orthogonality conditions and we ignore them in our discussion 8 . 

The most important fact is that the constraints can be explicitly solved. First, we write 
the solutions of the Liouville equations 1)21(1 in the form suggested by the conformal symmetry 
properties of the Liouville equation [TB"] . 



-W2 = a n {u)b n (v) - -g n a n (u)b n (v) = X n (u,v) 



(28) 



where the chiral fields a n (u), b n (v), a n {u) and b n (v) satisfy the equations (do not mix a n (u) 
with a n used above): 



a n (u)a' n (u) - a' n (u)a n {u) = 1, b n (v)b' n (v) - b' n (v)b n (v) = 1. 
Using (|29|) we can express a and b in terms of a and b and thus write X n as 



X n (u,v) = a n (u)b n (v) 



du 



du 



(29) 



(30) 



a 2 n (u) J bl(v)_ 

It is not so straightforward but not very difficult to rewrite the constraints (|2Uj) in the form: 



N 



a n (u) ^ b n {v) 

We first note that from (|18|) and (|2"3*)) it is easy to find the identities 

N N 
- diFdw + Vn) 2 = £ ln {di q n f . 

n=3 



(31) 



n=l 



It follows that 



N N 
C l = Y. q n ? - 28fq n ) = 4 £ 7n dfX n /X n 

n=l n=l 



Now, from the definition of X n it is not difficult to see that 

8 2 u X n <(n) dlX n &» 



Xn 



a n (u) 



Xn 



b n (v) 



(32) 



(33) 



(34) 



Using the relation In = we can explicitly solve these constraints. We show this for C u 
as the derivation for C v is similar. Let us introduce temporary notation a' n (u)/a n (u) = r n (u), 
use it in expression for C u and then shift all r n by an unknown auxiliary function R(u), i.e. 
define p n (u) = r n {u) + R(u). Thus we may express C u in terms of p n {u) and R(u): 



A' 



± C - V-v an " {u) 
4 n=l 7B «»(«) 



A" 



A 



(35) 



71=1 



71=1 



where we used the identity J2 In = 0. Now it is easy to see that the constraint C u = AD U 
will be solved if we take 



R(u) 



A 



Tn(pn(«) + Pl( u )) ~ D u{ 



u 



.n=l 



N 

^2jnPn(u) 
J Ln=l 



(36) 



In fact, the coupling constants have nothing to do with integrability and may be arbitrary numbers. 
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where p n (u) are arbitrary functions. Thus, if we choose a' n (u)/a n {u) = p n — R(u), where R 
is given by lj5BJl. the constraint C u = 4D U will be satisfied 9 . 
Now we introduce new moduli fields 



Pm(u) = p m {u) - - 



■ N 

E- 2 

n=l 



lnP 2 n {u) - D u (u) 



■ N 

E TnPr, 

.71=1 



(37) 



Using this definition and the identity j n = 0, it is easy to check that 

YllnVnW) = ^2 lnPn (u) , ^2lnP 2 n {u) = D u (u) = -^2{a k {u)f (38) 

(here and in what follows the limits of summation are omitted when the summation extends 
over all possible values of n = 1, ...,N and k = 1, ...,K). Repeating this derivation for the 
second constraint C v and defining the u-analogue, v n (v), of the moduli p, n (u), we finally write 
both constraints in terms of arbitrary chiral fields p n (u), v n (v), afc(u), fik{v): 

E 7«^(«) + \ EK W) 2 = o , E wl(y) + 1 E(^(^)) 2 = • (39) 

These constraints are equivalent to the original ones, (J20J), if the fields /, ip, ip, x are solutions 
of the equations of motion. 

Returning to the original moduli fields p n and to their relation to a' n (u)/a n (u) = r n {u) 
(and to their u-analogues) , one can see that the constraints are equivalent to the first order 
differential equations for a n {u) and b n {v): 

— — = p n (u)--— -— , —— = u n (v)--— -— , (40) 

a n {u) ^2^1nPn{u) b n {v) 2 LlnVniv) 

where p n (u) and v n (v) are arbitrary functions satisfying the constraints (J39|). This is a small 
'miracle' allowing us to explicitly solve the equations of motion (including the constraints) in 
terms of a 'sufficient' number of arbitrary functions. The number of independent arbitrary 
chiral functions in r.h.s. of (|4U|) is 2(A^ — 1 + K), where N — 2 is the number of the scalar 
matter fields ip, and K is the number of the free scalar matter fields \- Together with 
the 'gravitational' degrees of freedom 93, F (or ^i,^) the theory without constraints would 
have 2(iV + K) chiral degrees of freedom. We will show in a moment that the residual 
coordinate transformations allow us to further reduce the number of independent arbitrary 
chiral functions to 2(JV - 2 + K). 

By integrating the first order differential equations (|4*U|) for a n (u) and b n (v) we find the 
general solution of the A^-Liouville dilaton gravity in terms of the chiral moduli fields p n (u) 
and u n {v) satisfying the constraints ()39|) : 

a n (u) = lE^mMml 2 exp J dup n (u), b n (v) = |^7 m f m | 2 exp J dvv n (v). (41) 

Let us choose arbitrary chiral x-fields a^{n) and (3k( v )- Then the moduli fields p n (u) and 
v n {v) are not independent due to the constraints (|39]> . In addition, we may use residual coor- 
dinate transformation u — > U{u) and v —*V(v) to choose two gauge (coordinate) conditions. 
We will show in a moment that writing 

lnPn{u) |= U'(u) , I Y.Wn(v) \= V'(v) (42) 



E 



9 Note that the p„ — R are not independent functions. In fact, the chiral fields a n depend on iV — 1 arbitrary, 
independent functions of u. This will be clear in a moment from a somewhat different representation of the 
solution of the constraint. 
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is indeed equivalent to choosing (U,V) as a new coordinate system. With this aim, define 
new chiral fields, 

A n = exp J dufj, n (u) , B n = exp / dvv n (v) , (43) 

that we consider as the functions of the new coordinates (U, V) (due to implicit transforma- 
tions (g2J)). Using (jUJ) and (@2J) we may rewrite (J3UJ) in terms of A n (U), B n (V): 



Y n (U,V) = A n (U)B n {V) 
where we have defined 



1 r dU f dV 



Al(U) J B*(V) 



(44) 



Y n (U,V) = [U'(u)V'(v)}2X n (u,v) (45) 
Now, with the above definitions we may derive the metric /, 

N 

f(u, v) = U'{u)V'{v) J] [Y n (U, V)]^ a - = f(U, V)U'(u)V'(v), (46) 

n=l 

the dilaton (p, and the scalar fields ip m (m > 3): 

N N 

e^= II i Y n(U, V)}^\ = J] [Y n (U, v)]~ 2amnln . (47) 

n=l n=l 

We see that ds 2 = —4f(u,v)dudv = —4f(U,V)dUdV and thus everything is expressed in 
terms of the new coordinates (U, V). 

The constraints ()39|) for the moduli parameters can easily be solved by expressing one 
of the moduli in terms of others. However, it may be more convenient and instructive to 
introduce new moduli that make clear the topological nature of the moduli space and thus 
may allow interesting topological classification of the solutions of the iV-Liouville theory. To 
simplify notation we define the new moduli for the case of zero x-helds (the general case is 
not essentially different). Thus, setting x = 0we introduce the following unit (iV — l)-vectors 
(recall that 71 < and 7^ > for k > 2): 

&(«) = — r , rj k (v) = — r , k = 2, N . 48) 

These vectors, moving on the (N — 2)-dimensional unit sphere S^ N ~ 2 ^ , determine the solution 
up to a choice of the coordinate system, which can be fixed by the above gauge conditions 
(|42|) , They now look as follows: 

U'(u) = | W («)|(1- cos V'(v) = \ W («)|(1- cos 6 v (v)), (49) 

N N 

cos 9^(u) = %ik(u) = 7I, cos9 v (v) = ^ = IV, (50) 

k=2 k=2 

where 7 is the constant unit (N — l)-vector, 7^ = (t^/Iti 1)^ • 

The moduli fields £(u) and f]{u) are very useful because they give very simple representa- 
tion of the most important solutions and, in particular, visualize relations between solutions 
of different dimensions. For example, if the vectors £(u) and r)(v) are constant and equal 
(£ = v)i they give a static or a cosmological solution. The static solution has a horizon if 
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£ = fj = 7 (see PQ). Another interesting static solution, which is flat at infinity, is given 
by the constant vectors satisfying the condition £ = i) = 5, where 5k = (7fc a fe/lTi l a i)^ (it 
is not difficult to check that this vector is indeed of unit length). Two - dimensional solu- 
tions may be represented by pairs of curves in S^ N ~ 2 ^ that we denote by (£(m), v( v ))- Those 
that asymptotically interpolate between the above one - dimensional solutions are of special 
importance and define additional structure on S( N ~ 2 \ which may be used for a physically 
motivated classification of the solutions of the iV-Liouville theory. 

Especially interesting are the new solutions that, probably, may be localized in space 
('lumps' or soliton - like waves of scalar matter). They correspond to constant but not equal 
vectors ^ and f) (obviously, this is possible for N > 3). In the coordinates (U, V) (U + V = r, 
U — V = t) the solutions may be written as 

Y n (U,V) = C nC osYi[R- 1 ({r-r n )-v n t)], (51) 

where the parameters C n , R n ,r n ,v n rationally depend od the moduli and other parameters 
characterizing the Lagrangian and the coordinate system (recall that this parameters are not 
independent, due to the constraints and the gauge conditions). Although the exponentials 
e q = Y~ 2 are localized in space, this is not necessarily sufficient for the matter fields to 
be localized. However, using freedom in choosing moduli and other parameters it may be 
possible to write a truly localized solution. A detailed derivation and investigation of these 
interesting solutions will be presented in a separate publication. 

We thus have the general solution of the 1+1 dimensional dilaton gravity coupled to 
any number of scalar fields. It is explicitly expressed in terms of a sufficient number of 
arbitrary chiral fields 10 and thus we may solve the Cauchy problem and study the evolution 
of cosmological or black hole type solutions, etc. The representation of the general solution 
in terms of the chiral fields a n {u) and b n (v) may give us a good starting point in attempts 
to quantize our A^-Liouville dilaton gravity. Even more useful may be the chiral moduli 
fields n n (u) and v n (v) (or and r)k{v)). In terms of these moduli fields the dimensional 

reduction of the solutions becomes very transparent and this may simplify the derivation 
and physical interpretation of the evolution of one - dimensional solutions and suggest new 
approaches to quantization based on analogy with the simple 1-dimensional case 11 . 

4 Discussion and outlook 

The explicitly analytically integrable models presented here may be of interest for different 
applications. Most obviously we may use them to construct first approximations to generally 
non integrable theories. Realistic theories describing black holes and cosmologies are usually 
not integrable. However, explicit general solutions of the integrable approximations may 
allow one to construct different sorts of perturbation theories. 

10 At first sight, the expressions for the fields /, <p, ip m may look nonlocal. However, denoting the indefinite 
integrals in 1441 by new chiral fields, we see that Y n are essentially local functions of new fields and of their 
derivatives. 

llr The one - dimensional models obtained by simple dimensional reduction may easily be quantized. Our 
formulation is as close to the one - dimensional theory as possible and this supports hopes to finally find a 
quantum version of the iV-Liouville theory. It should be emphasized that, up to now, the interesting results 
obtained in quantum Liouville theory did not shed much light on the iV-Liouville theory. This may be explained 
by the fact that the main difficulty and the main content of the iV-Liouville, as distinct from the standard 
Liouville theory, lies in the constraints. 
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For example, spherically symmetric static black holes non minimally coupled to scalar 
fields are described by the integrable 0+1 dimensional iV-Liouville model. Similarly, spher- 
ically symmetric cosmological models may be described by 1+0 dimensional 7V-Liouville 
theory. However, the corresponding 1+1 theory is not integrable because the scalar coupling 
potentials Z n are not constant (actually, Z n ~ ip). To obtain approximate analytic solutions 
of the 1+1 theory one then may try to approximate Z n by properly chosen constants. 

This approach may be combined with the recently proposed analytic perturbation theory 
allowing one to find solutions close to horizons for the most general non integrable 0+1 dilaton 
gravity theories |16j . Near the horizons we can use the integrable 1+1 dimensional dilaton 
gravity (with Z n = — 1) as a good approximation to a realistic theory (with Z n depending on 
the dilaton p). 

In cosmological applications, the behaviour of the 1+1 dimensional solutions near the 
singularity at <p = Ois of great interest (see e.g. |17|). Integrable 1+1 dimensional iV-Liouville 
theories could give, at best, a rough qualitative approximation of the exact solutions near the 
singularity. A more quantitative approximation might be obtained by first asymptotically 
solving the exact theory in the vicinity of p = and then sewing the asymptotic solutions 
with those of the integrable theory. To realize such a program one needs a very simple and 
explicit analytic solutions of the integrable theory. Our simple model having the solutions 
represented in terms of the moduli £ and f\ may give a good starting point for such a work. 
Of course, before applications to realistic cosmologies become possible, one should study in 
detail and completely classify and interpret the behaviour of the 1+1 dimensional solutions 
and their precise relation to the 1+0 dimensional reduction. 

The reduction from dimension 1+1 both to dimension 1+0 ('cosmological') and to dimen- 
sion 0+1 ('static' or 'black hole') is especially transparent in the moduli representation for the 
solutions of the 1+1 dimensional iV-Liouville model. However, as emphasized in [T], the whole 
procedure of the dimensional reduction should be reconsidered from a more general point of 
view, taking into account more general dimensional reductions. A detailed consideration of 
generalized dimensional reductions will be given elsewhere. 

5 Appendix 

Here we sketch a simple approach to solving the pseudo orthogonality conditions for a mn 
with additional restrictions a\ n = 1 + a n , am = ^ — a n . The equations one has to solve are 



These N(N — l)/2 equations are nonlinear but there exist a simple recursive algorithm re- 
ducing their solution to solving linear equations. To find it one has to choose which of the 
parameters should be regarded as unknown ones. Analyzing the cases N = 3 and N = 4 
one may find that the convenient choice of the division of the parameters into unknown and 
arbitrarily fixed is the following. Let us fix 033 and a mn for m > n, m > 3 and thus all other 
parameters are unknown. First, consider the equations for i,j = 1,2,3. Denoting the r.h.s. 
of ((52|) by Ak where (ijk) = (123) cyc ii c we find 




n=3 



(52) 



2ai = -At + Aj + A k , (ijk) = (123) cyc i ic . 



(53) 
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The next step is to consider the equations (|52*|) for i = 1,2,3 and j = 4. These constitute 
three linear equations for three unknowns 04, 034, 044. They have a unique solution provided 
that the r.h.s. and the determinant A do not vanish. 

Now the general construction should be clear. On the j-th step, where j = 4, ...,iV, we 
have j — 1 equations for j — 1 unknowns, a,j and a n j for 3 < n < j: 

3 N 

- 2aj + ^2 a ni a nj = 2a,i - ^ a ni a nj , i = l,...,j-l. (54) 

n=3 n=j+l 

The r.h.s. of these equations depend both on arbitrary and previously found parameters. 
The nontrivial part of the procedure is that the determinant, starting from j = 5 depend on 
the previously found parameters and thus the condition A 7^ is difficult to control. For 
j = 4 the determinant, 

A = 2[(a 33 - a 3 i)(a 43 - a 42 ) - (a 33 - a 32 )(a 43 - a 4 i)] , 

depends on the arbitrary parameters only. 

Finally, let us write the parameters of the general model with N = 3: 

2m = a 3i a 3j - a 3j a 3k + a 3k a 3i , 7," 1 = (a 3i - a 3j )(a 3i - a 3k ) , (ijk) = (123) cyc i ic . 
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